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ABSTRACT 

Despite the efforts of the past decade, the origin of the bimodal horizontal- 
branch (HB) found in some globular clusters (GCs) remains a conundrum. In- 
spired by the discovery of multiple stellar populations in the most massive Galac- 
tic GC, uj Centauri, we investigate the possibility that two distinct populations 
may coexist and are responsible for the bimodal HBs in the third and fifth bright- 
est GCs, NGC 6388 and NGC 6441. Using the population synthesis technique, 
we examine two different chemical "self-enrichment" hypotheses in which a pri- 
mordial GC was sufficiently massive to contain two or more distinct populations 
as suggested by the populations found in u Gen: (1) the age-metallicity relation 
scenario in which two populations with different metallicity and age coexist, fol- 
lowing an internal age-metallicity relation, and (2) the super-helium-rich scenario 
in which GCs contain a certain fraction of helium-enhanced stars, for instance, 
the second generation stars formed from the helium-enriched ejecta of the first. 
The comparative study indicates that the detailed color-magnitude diagram mor- 
phologies and the properties of the RR Lyrae variables in NGC 6388 and NGC 
6441 support the latter scenario; i.e., the model which assumes a minor fraction 
(~ 15 %) of helium-excess (Y ~ 0.3) stars. The results suggest that helium 
content is the main driver behind the HB bimodality found most often in mas- 
sive GCs. If confirmed, the GC-to-GC variation of helium abundance should be 
considered a local effect, further supporting the argument that age is the global 
second parameter of HB morphology. 
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Subject headings: Galaxy: formation — globular clusters: general — globular 
clusters: individual (NGC 6388, NGC 6441) — stars: horizontal-branch — stars: 
RR Lyrae variables 



INTRODUCTION 



The horizontal-branch (HB) morphology of a globular cluster (GC) is defined as the color 
distribution of its HB stars. The physical cause of the wide diversity in the HB morphology 
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though several other variables may simultaneously play roles in H B morphology, the idea 
that GC age is the major second parame ter has been very popular (ILee. Demarque. &: Zinn 
I994J : iRey et al.ll200ll : lYoon fc Leell2002l ). However, one of the strongest arguments against 
the age hypothesis ha s been the existence of a number of G Cs whose HB color distribution 
is distinctly bimodal (IRood et al.lll993l : IStetson et al.lll996l ). An explanation based purely 
upon age as the second parameter suggests that, at a given metaUicity, the population as- 
sociated with the red HB is a few gigayears younger than that associated with the blue HB. 
This cannot be reconciled with the traditional "single-population" picture of GCs. 

Perhaps the most striking examples of HB bimo dality are the two Galactic bulge GCs , 
NGC 6388 and NGC 6441, along with NGC 2808 JPiotto et al.l I2OO2I : ISosin et al.l E997l ). 
NGC 6388 and NGC 6441 are relativ ely metal-rich GCs with [Fe/H] = -0.60 and -0.53, re- 
spectively (lArmandroff fc Zinrull998l ). Observations from the Hubble Space Tele scope (HST) 



have revealed the pres e nce of a significant fraction of extended blue-HB stars ([Piotto et ah 



Pritzl et al. 


2003: 


Piotto et al. 


2002 



2004 



Catelan et al.ll2006h . in addition to a majority of red-HB stars as would be expected given 



their metallicity. 

Recent studies suggest that the long-standing puzzle of the HB bimodality may no 
longer be a complete mystery. It has been discovere d that at least four discrete populations 



2006 



coexist in the most massive Galactic GC, u Gen (iLee et al.l Il999bl : iStandford et al 
and references therein) . This finding obviously opposes the conventional "single-population" 
picture of GCs, and also provides an instructive precedent for the bimodal-HB feature in 
other GCs. The multiple populations in u Gen imply an internal age-met allicity relation 
(AMR), in that stars having a higher metallicity are younger. Furthermore, subsequent 
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observations have revealed the presence of a marked double main-sequence (MS) with a 
minority popula t ion of bluer and /o r f ainter MS sta rs, w hich are separated fr om the majority 
flAndersonlEo02l : ISedin et al.ll2004h . iNorrisI tooi ) and IPiotto et al.l J2005h have suggested 
that a greatly enhanced helium abundance (AY = 0.12 ~ 0.14) can explain the split MS of 
LJ Cen. Lee et al. (2005) have further shown that the helium-enhanced population in the 
MS is most likely the progenitor of the extreme HB (EHB) stars. 

It has also corne to our attention that the integrated luminosities of NGC 6388 and 
NGC 6441 (iHarrid 119961 ) make them th e third and fifth brighte s t GCs next to u Cen and 
M54 among the ~ 150 Galactic GCs JYoou. Lee, fc Leel I2OO0I : iRee et al.l [2002! ). We are 
particularly inspired by the fact that five of the 10 brightest GCs can be characterized by 
the "composite" nature of color-magnit ude diagrams (CMDs ), such as multiple red-giant- 
branches (RGBs) and/or biniodal HBs (Ivan den Berglj|l996l ). The five GCs in this group 
are u Cen, M54 (ISiegel et al.l 120071 ). the two GCs of interest in this study, NGC 6388 and 
NGC 6441, and the sev enth br i ghtest GC, NGC 2808 (Table 1, 6th column). On the other 
hand, according to the iHarrid (119961 ) catalog, the 11 confirmed or suspected bimodal-HB 
GCs ( ICatelan et al.lll998l ) are always brighter than the turnover of the l^-band luminosity 
function {My* ~ —7) of the Galactic GC system (Fig. la). The relative fraction of the 
combined sample of the bimodal-HB GCs, u Cen, and M54 progressively increases at the 
higher brightness bin (Fig. lb). It may be that there is a link between the massive bimodal- 
HB GCs and the most massive multiple-population GCs - u Cen and M54. 

The purpose of this paper is to explore, with the stellar population synthesis technique, 
the possibility that two distinct stellar populations may coexist and be responsible for the 
HB bimodality in NGC 6388 and NGC 6441. We assume the two chemical "self-enrichment" 
processes that were proposed to explain discrete subpopulations in u Cen: (1) the two 
hypothetical subpopulations in NGC 6388 and NGC 6441 follow an internal AMR (hereafter 
"the AMR scenario"), and (2) NGC 6388 and NGC 6441 contain a certain fraction of super- 
helium- rich (SHR) stars (hereafter "the SHR scenario"). Further investigation is in progress 
to determine whether one of these scenarios can account for other GCs thought to contain 
bimodal HBs JCatelan et al.lll998h . including NGC 1851, NGC 2808, and NGC 6229. 



The paper is organized as follows. A brief description of the stellar population models 
is given in §2. We compare the synthetic CMDs (§3) and the model RR Lyraes (§4) to 
observed data. We argue that the bimodality of HBs and the unusual properties of the 
RR Lyrae variables observed in NGC 6388 and NGC 6441 are better explained by the SHR 
scenario. Finally, we discuss the implications of our results in §5. 
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STELLAR POPULATION MODELS 



The present models are constructed using the Yonsei Evolutionary Population Synthe- 



sis (YEPS) code ( 


i'ark & Lee 


1997; 


Lee et al. 


1999b: 


Lee. Yoon & Lee 


2000; 


Rey et al. 


2001; 


Yoon & Lee 


2002 


Lee et al. 


2OO5I; 


Yoon. Yi. & Lee 


2006 


). The models use a set of stellar 



libraries for the MS-RGB and post-RGB evolutionary tracks that were built using the same 
input physics for consistency. Combined with RR Lyrae pulsation th eories, the synthetic pop- 
ulation models can pred ict the properties of the RR Lyrae variables (ILee. Demarque. fc Zinn 
I990I : lYoon fc Leell2002l ). It is not our intention to discuss the population synthesis technique 
in this paper, and readers are referred to the above papers and Yoon et al. (2008, in prep.) 
for the details of the YEPS model. Table 2 and Table 3 summarize the model's ingredients 
and input parameters employed in this study. 



3. COLOR-MAGNITUDE DIAGRAMS 



3.1. The Age-Metallicity Relation (AMR) Model 

Figure 2 co mpares our popu lation models to the observed HST CMDs of NGC 6388 
and NGC 6441 (Rich et al. 1997). We note that more recent HST CMDs now exist for the 



GCs flPritzl et al 



2003 



Catelan et al. 2006). For NGC 6441, however, we selected Rich et 



al.'s CMD instead of that of iPritzl et al.l (120031 ) as the former is better fit for the purposes 
of our CMD analysis in this section. The snapshot study by Pritzl et al. is abundant in 
information on variable stars and clearly more appropriate for the RR Lyrae analysis, so we 
fully exploit their results in §4. For the sake of consistency ^ yith the case o f NGC 6441, we 
also opt for Rich et al.'s CMD over the superseding CMD by lCatelan et al.l (120061 ) for NGC 
6388. 

To simulate CMDs based on the AMR hypothesis (Figs. 2b & 2e), we first generate 
model CMDs assuming a common age of 13 Gyr, but different [Fe/H] of —0.60 for NGC 6388 
and —0.53 for NGC 6441. Once the model shows reasonable agreement with the MS through 
the RGB and red portions of the HB, we then, according to our working hypothesis, add the 
second minor component following a mild internal AMR. We adjust the AMR parameters 
until the model CMDs mimic both the overall appearance and the population number ratios 
of the HB. Given the level of photometric accuracy, the purpose of our CMD fit is neither to 
determine the absolute age of the GCs nor to make a definitive statement about their exact 
AMR parameters. For a direct comparison with observations, we carried out Monte Carlo 
error simulations based on the actual observational uncertainties. 
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The occurrence of a rather strong differen tial reddening across the fields of NGC 6388 
and NGC 6441 is a well-known phenomenon ((piotto et al. 1997 : Heitsch fc Richtler 1999 



Layden et al.lll999l : iRaimondo et al.l l2002l : iLaw et al.ll2003l ). Interestingly, the differential 
reddening effect facilitates the development of markedly sloped red clumps. By including 
such an effect, our model can reproduce the tilted red clump found in the GCs. However, 
several studies have argued that the observed differential red dening would not be enough to 
turn a "normal" red HB into a significantly slope d structure (ICatelan fc de Freitas Pacheco 
19961 : iBeaulieu et al.l 1200 ll : IRaimondo et al.ll2002l ). It is thus important to assess the extent 
to which the models are able to reproduce the red clumps. 

Figure 3 illustrates our numerical experiments on the impact of differential reddening 
upon these structures. The middle panels (Figs. 3b & 3e) show the flat red clumps without 
differential reddening, and the right panels (Figs. 3c & 3/) show the tilt of the clumps due 
to this effect. The simulations suggest that a differential reddening value of aE{B-v) = 0.03 
(mag.) is sufficient for the models to reproduce the observed tilt of the red clumps. This is 



compa rable to the observed value of AE{B — V^) < 0.1 (~ 3 x c^g^R_^/^ l in IRaimondo et al. 
(l2002h. and slightlv small e r than AE(B-V) 0.15 f~ 5 x aF.(p -v\) in [PiottoelalJ fll997l ). 



Heitsch fc Richtled (119991 ) , iLayden et al.l (119991 ) , and iLaw et al.l (120031 ) . Figure 3 also shows 
that the models reproduce the observed RGB bump positions. It should be stressed that, 
in addition to the shape of red HB clumps, the differential reddening effect seems to be 
necessary to reproduce the observed scatter and possible tilt of the RGB bumps. 

We now turn back to Figure 2. The comparison of the CMDs between observed data 
(Figs. 2a & 2d) and the models (Figs. 2b & 2e) shows reasonable matches from the MS 
through the HB. The colors of RGB stars are sensitive to metallicity, especially in the high- 
[Fe/H] regime. Yet, the metallicity differences between the two hypothetical populations 
are small enough to produce apparently single RGBs within observational uncertainties. 
Decreasing age shifts the color of the RGB towards blue, and this effect, albeit small, also 
helps to generate the single RGBs. 

Another important result illustrated in Figure 2 is that the observed bimodal HBs 
can be reproduced by the sum of the red HB (from the metal-rich younger component) 
and the blue HB (from the metal-poor older component). Interestingly, theories of stellar 
evolution predict that increasing age or decreasing metallici ty yields a bluer HB, wherea s 
decreasing age or increasing metallicity produces a redder HB (ILee. Demarque. fc Zinnlll994l ). 
Hence, it is obvious that an internal AMR splits readily t he HB of each population in two 



leading to bimoda l HBs. Moreover, our population models (ILee et al.lll999al : iRey et al.ll2001 



Yoon &: Ledl2002l ) show that the HB mor phology is more sensitive to GC age by a factor 
of two when compared to earlier studies (ILee. Demarque. &: ZinnI 119941 ) . This results in a 
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substantial reduction in the required age difference between blue- and red-HB populations . 
This is significant because all of the bimodal-HB GCs found to d ate (ICatelan et al.lll998l ) 
have no noticeable structure in the CMDs near the MS turn-offs (jPiotto et al.ll2002l ). thus 
leaving little room for an age spread of more than a couple of gigayears. The difference in 
[Fe/H] and the ages required to reproduce their CMD morphologies are 0.35 dex and 2.0 
Gyr for NGC 6388, and 0.22 dex and 2.0 Gyr for NGC 6441 (Table 3). 

There are several versions of the semi-empirical mass-loss formula used for reproducing 
the obse rved amount of themass loss in giant stars. Various mass-loss formulae are summa- 
rized by ICatelanI J2000l . boosh . We note that a different choice of formulae somewhat alters 
the required difference in [Fe/H] and age between the b lue and r e d HB s, bu t does no t affect 
our conclusion. For instance, according to Figure 11 of ICatelanI (120001 ). the iReimerd (Il975l ) 
formula would allow Aage = 2.0 Gyr centered at 13 Gyr to yield a 5 % difference in AM 



(i.e., t he mass loss along the RGB). If one adopts Equation A5 (the steepest one) in lCatelan 
( I2OOOI ). then 5 % in AM translates into Aage = 2.3 Gyr centered at 13 Gyr. Thus, even the 
most extreme choice of the mass-loss formula would cause only a ~ 15 % [= (2.3— 2.0)/2.0] 
change in the age difference between red- and blue-HB populations. In order to ins pect the 
metallicity dependence of the formulae, we take the steepest example in Figure 4 of ICatelan 
(120051 ) instead of the Reimers' formula, and find that the required difference in [Fe/H] be- 
tween red- and blue-HB populations in the AMR model becomes smaller by about a factor 
of eight. This makes the reproduction of bimodal HBs even easier. 

The above analyses lead us to conclude that a mild internal AMR between the two 
hypothetical populations may be a plausible cause of the HB bimodality, together with the 
single MSs and RGBs within observational uncertainties. 



3.2. The Super-Helium-Rich (SHR) Model 

We now consider the SHR hypothesis. The model CMDs in Figures 2c & 2/ are con- 
structed based on this assumption. To simulate the CMDs, we adopt [Fe/H] = —0.60 for 
NGC 6388 and —0.53 for NGC 6441, with a common age of 13 Gyr and a helium abundance 
of Y = 0.245. The major dominant populations are the same in both the AMR (the metal- 
rich younger component denoted by "1" in Figs. 26 & 2e ) and the SHR (the normal-helium 
component denoted again by "1" in Fig. 2c & 2/) models. We then add a helium-enhanced 
population according to our working hypothesis. The helium abundance and the number of 
stars are adjusted until the best match between the modeled and observed CMDs is achieved, 
especially for the blue HB appearance. Observational errors and differential reddening effects 
are also simulated in the same manner as in the AMR model. 
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Piotto et al.l (120051 ) report that the stars of the bluer component of the spht MS in uj Cen 
are more metal-rich (by 0.3 dex in [Fe/H]) than those of the dominant redder component. 
For NGC 6388 and NGC 6441, however, there is no such constraint dictated by observations. 
For a reference, in order for a 0.3 dex more metal-rich population to mimic the observed 
blue-HB structure in NGC 6388 and NGC 6441, the required increase in helium abundance 
(AY) is as small as 0.02. Hence, in the SHR scenario, we consider the helium-enriched 
population to have the same [Fe/H] as the normal-Y population. The parameters required 
to reproduce their CMD morphologies are listed in Table 3. 

The model CMDs based on the SHR assumption match well with observations from the 
MS through the HB. The variation in helium abundance has a relatively weak effect on the 
MS to the RGB, so the simulations show no indication of bifurcation on the MS through 



the R GB in each mod el CMD. The mul tiple MS feature, as found in uj Cen ( iPiotto et al. 



20051 ) and NGC 2808 (IPiotto et al.l 120071 ). is not feasible in the available CMDs which con- 
tain considerable observational uncertainty under the MS turnoffs. More importantly, the 
observed bimodal HBs can be reproduced by the sum of the red HB (from the normal helium 
component) and the blue HB (from the helium-enriched component). For identical values 
of the total mass a nd helium core rnass, theoretical models show that the HB stars with a 
higher Y are bluer jRood et al.lll97ol : ISweigart fc Grosslll976h . Moreover, for a given age and 
metal abundance, a helium-enriched HB star has a thinner hydrogen envelope surrounding 
the helium-burning core, making the HB even bluer. As a result, bimodal HBs with single 
RGBs can be achieved given the assumption that ~ 15 % of the stars in NGC 6388 and 
NGC 6441 have an enhanced helium abundance of AY ~ 0.05 (Table 3). 

The CMD morphology analysis suggests that both the AMR and SHR models are gen- 
erally in first-order agreement with observations, reproducing the bimodal HBs along with 
the single RGBs. We proceed to make a detailed comparison of the HB structures in the 
two models. 



3.3. Comparison between the AMR and SHR Models 



Apart from the overall agreement in the CMD morphology between the observations and 
the models, there is an important feature of HBs that must be explained in detail. The ob- 
served HBs seem to slope upward with decreasing {B — V), and their upper parts are brighter 



than the bulk of the red HBs by ^ 0.5 inag in V ( Rich et al 



Raimondo et al.ll2002l : IPritzl et al.ll2003l : ICatelan et al.ll2006f ). 



19971 : ISweigart fc Catelanlll998 



The AMR models in Figures 26 & 2e can only partially recreate the feature by virtue of 



- 8 - 



the HB evolutionary effect (iLee. Demarque. fc Zinrull990l ). If there are two HB populations 
as predicted by the AMR model, the stars at the top of the blue HB are most likely to be 
highly evolved and thus brighter than stars near the zero-age HB (ZAHB). As a result, the 
HB evolutionary effect promotes the difference in brightness between the upper part of the 
blue HB and the red clump. However, inspection of Figures 2a & 2d suggests that the RR 
Lyraes in NGC 6388 and NGC 6441 appear to be dominated by stars near the ZAHB. It is 
therefore unlikely that the brightness of the blue HBs of NGC 6388 and NGC 6441 support 
the AMR scenario. 

In contrast to the AMR model, the interesting HB appearance can be more readily 
reproduced by the SHR model (Figs. 2c & 2/), because helium-rich HB stars are intrinsically 
bright er. In an attempt to expla in the sloping blue HB extensions in NGC 6388 and NGC 



6441, ISweigart fc Gate 



However, 



Layden et al. 



an! (I19981) have first proposed a high helium abundance (Y ~ 0.4). 



(|l999|) claimed that the estimated Y values of NGC 6388 and NGC 
6441 via the R-method strongly disfavor that scenario. Layden et al. showed that the models 
in which Y = 0.38 and 0.43, presented in Sweigart & Catelan, inevitably yield R = 3.4 and 
3.9. These are about 3 a away from the observed value of R = 1.6 ± 0.7, which corresponds 



to Y 



2S+° °^ 
"-"•^^-o.os- 



Although our SHR model in this study assumes a subpopulation with a 
high helium abundance (Y ^ 0.3), its number fraction is as low as ~ 15 %. A calculation 
shows that the number-weighted mean Y of the GCs in our SHR model is as low as Y ~ 
0.254 for NGC 648 8 and 0.252 for NGC 6441, in complete accordance with 0.25l2;[]^ by 
Lavden et all Jl999h . 



Observations indicate that NGC 6388 and NGC 6441 contain a small portion of ad- 
ditional EHB stars on their HBs. For the AMR model to reproduce this detail, it should 
employ a unreasonably large mass dispersion (ctm ~ 0.06 Mq) on the blue HBs. This is three 
times larger than the commonly used v alue of aM = 0.02 Mq ( iLee. Demarque. fc Zinnlll994j : 



Catelan et al.lll998l : lYoon fc Ledl2002l ). which is currently used for the red HBs in our mod- 
els. In contrast, the SHR model requires only ctm — 0.03 Mq to simulate the EHB structure. 
This is because, for larger helium abundances, the blue loops in the evolutioiiary tr acks can 
become considerably longer, reaching higher effective temperatures (jSweigartl 119871 ). 



Alternatively, there is a possibility of the presence of a minor third subpopulation which 
is responsible for the EHBs in NGC 6388 and NGC 6441. In the AMR hypothesis, the third 
population presumably has an older age than the underlying metal-poor old population. Our 
model suggests an age as high as > 17 Gyr. Such an unrealistically old age would not be 
acceptable. In the SHR model, on the other hand, the third population should have a larger 
AY than the blue-HB component. Our present model gives a helium abundance of Y ~ 0.33 
- 0.34 for the EHB stars in NGC 6388 and NGC 6441. 
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It is interesting to note that the SHR scenario is more consistent with EHB (a. k. 
"blue hook") stars that are known to be fainter than the redder HB stars and below 



the ZAHB locus. 



Lee et all f l2005f ) show that the hottest EHB stars found in NGC 2808 



( iBrown et al.ll200ll ) are fainter th a n the redder HB stars as a natural consequence of their 
high Y, consistent with ISweigartI (119871 ) predictions. This may indicate that a high Y is 
sufficient per se to explain the hottest blue- hook sta rs without necessarily inv oking a delayed 
helium flash (ILanz et al.ll2004j : iMoehler et al.ll2004j ). Readers are referred to ICatelanI (120051 ) 
for various characteristics of EHB stars. The origin of EHB stars and their possible link 
to the HB bimodality are interesting issues worth further investigation. In this study we 
consider the EHB stars as blue HB stars when calculating the number ratios of blue and red 
HB stars. 



Clementini et al.l (120051 ) have shown that the mean metallicity of the RR Lyrae stars in 
NGC 6441 is c lose e nough to the typical metal abundance for this cluster. More recently, 
Gratton et al.l (120071 ) have found no clear sign of star-to-star Fe abundance variation in NGC 
6441 from the Giraffe spectrograph at VLT2. Hence, there appears to be a growing body 
of evidence that contradicts t he notion that RR L yraes and blue HB stars have a lower 
metal abundance. Moreover, iMoehler fc SweigartI (120061 ) have reported that the physical 
parameters of the cool blue HB stars in NGC 6388 are consistent with the predictions of the 
helium enrichment scenario, adding support to the SHR hypothesis. 

We have compared the AMR and SHR models in terms of the CMD morphology, and 
concluded that the SHR model is more successful at reproducing the details of the observed 
CMDs of NGC 6388 and NGC 6441. There is yet another important feature to be accounted 
for - their unusual RR Lyrae properties. This is the subject of the following section. 



RR LYRAE VARIABLE STARS 



Obseryatiori s ( Sweigart &: Catelanlll998 : Layden et al. 1999; Clement et al. 2001; Pritzl et al 



2000, 2001, 2002, 2003 



Corwin et al.ll2006l ) have revealed unusual properties of the RR Lyrae 



variables in NGC 6388 and NGC 6441. First, the mean pulsation periods of both the ah- and 
c-type variables ar e too long for the i r metallicities, so t hey do not fall into either of the usual 
Oosterhoff groups (IOosterhofBll939l : lYoon fc Leell2002[ ). Second, the number fraction of type 
c variables, N{c)/N{ab + c), falls between the two Oosterhoff groups, but closer to that of the 
group II than of the group I. This is unusual because the traditional group II GCs are always 
metal-poor ([Fe/H] < —1.6). These findings have posed yet another serious challenge to our 
understanding of the GC stellar population. It is interesting to see whether the present 
models based on the AMR and the SHR hypotheses reproduce the peculiar properties of RR 
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Lyrae variables in NGC 6388 and NGC 6441. 



4.1. The Pulsating Periods 



Figure 4 presents a direct comparison between the observed and predicted RR Lyraes 
in NGC 6388 and NGC 6441 in terms of their pulsating perio ds as a function of {B — V)q. In 
Figures 4a & 4d, we s how 12 RR Lyraes foun d in NGC 6388 JPritzl et aljl2002h . and 63 RR 
Lyraes in NGC 6441 JPritzl et al.lboOll . l2003l ). To obtain the fundamental periods (Pp), the 
c-type period (Pr) has been fundament alized using the eq uation, Log Pp = Log Pc + 0.13 
(ICastellani fc Quartal 119871 : iLee. Demarque. &: Zimj|l990l ). As shown in Figure 2, the RR 
Lyrae stars are predicted to belong to the blue-HB components in both the AMR and SHR 
models. Since the inferred physical parameters of the blue-HB populations of NGC 6388 and 
NGC 6441 are almost identical (Table 3), one may consider the RR Lyrae populations of the 
two GCs to be twins. Indeed, Figure 4a reveals that, despite a difference in the observed 
number of RR Lyraes in NGC 6388 and NGC 6441, their RR Lyrae distributions in the Pp 
vs. {B — V)o diagram are statistically indistinguishable from each other. On these grounds, 
we combine the catalogs of NGC 6388 (12 RR Lyraes) and NGC 6441 (63 RR Lyraes) in 
order to minimize the small number statistics. 

The AMR model predicts shorter periods at a given {B — V)o than observations show. 
Figures Ab & 4e show examples of the Monte Carlo simulations. The 10,000 Monte Carlo 
realizations under the AMR assumption give (Pf) = 0.555 ± 0.041 (day) from 12 ± 3 model 
RR Lyraes (NGC 6388) and 0.567 ± 0.017 (day) from 63 ± 7 (NGC 6441). Each of the (Pp) 
values is 2.6 a and 5.5 cr from the observed value of (-Pf) = 0.661 (day) for the combined 
NGC 6388 and NGC 6441 data. For the AMR model to achieve a longer mean period, more 
highly-evolved RR Lyraes and thus a bluer HB are required. However, the number of these 
evolved and bright RR Lyrae stars is expected to decrease sharply as the HB morphology gets 
bluer. This problem has been investigated in detail by iPritzl et al.l (120021 ) . who concluded 
that the model, which produces the observed ratio between blue HB stars and RR Lyraes, 
cannot reproduce the exceedingly long mean period of NGC 6388. Therefore, the AMR 
model appears to be in conflict with the available RR Lyrae observations. 

On the contrary, the SHR model appears to succeed in reproducing both the mean 
Pp and the observed {B — V)q-Pf correlations. Figures 4c & 4/ show examples of the 
Monte Carlo simulations. The 10,000 Monte Carlo realizations performed under the SHR 
assumption give (Pp) = 0.629 ± 0.041 (day) from 12 ± 3 model RR Lyraes (NGC 6388) 
and 0.645 ± 0.018 (day) from 63 ± 7 (NGC 6441). Each of the (Pp) values exhibits < 1 a 
agreement with the observed value of (Pp) = 0.661 (day). As in the AMR model, all of the 
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RR Lyraes belong to the blue-HB population in the SHR model. However, this model does 
not require the presence of highly evolved RR Lyraes to produce long-period RR Lyraes, as 
helium-enhanced RR Lyraes are intrinsically brighter and of longer periods. Since most of 
the RR Lyraes exist at the stage of the ZAHB, their significant numbers are also reproduced. 
We emphasize that AY = 0.05 reproduces simultaneously the peculiar RR Lyrae periods and 
the ratio between blue HE and RR Lyrae stars. 



Pritzl et al.l (120021 ) has attributed great importance to the presence of long-period c- 



type RR Lyraes in NGC 6388 and NGC 6441, which are rarely found in other GCs. The 
mean periods of c-type variables, {Pc), in NGC 6388 and NGC 6441 are 0.36 (day) and 0.38 
(day). Note that the average (Pc) for Galactic GCs is known to be 0.33 (day), respectively. 
In the AMR model, the assumed c-types (i.e., {B — V)o < 0.3) are concentrated toward the 
ZAHB (Figs. 46 & 4e) and thus have short periods. To include evolved c-types, the blue 
HB component must become bluer. However, the model shows tha t this shift g i ves ri se to 
a significant decrease in the number of RR Lyraes, as discussed by iPritzl et al.l (120021 ). In 
contrast to the AMR hypothesis, a raise in (Pc) agrees with the SHR scenario. The high 
helium content increases the zero-point of RR Lyrae brightness, leading to c-types of longer 
periods (Figs. 4c & 4/). 



4.2. The Oosterhoff Groups 



Figure 5 shows the correlation between the mean period of type ab RR Lyraes {{Pab)) 
i n GCs and thei r [Fe/H]. In Figure 5a, we display the well-known Oosterhoff dichotomy 
JOosterhofjE939l ). along with the unusual d i stribu tion of NGC 6388 a nd N GC 6441. The 
{Pab) data are obtained from I Clement et al.l (120011 ) and ICorwin et al.l (120061 ). One can see 
that the Galactic GCs are mainly divided into two distinct groups according to mean period 
and metal abundance: Oosterhoff group I (filled squares; (Pab) — 0.55 day, and [Fe/H] > 
-1.8) and group II (open squares; (Pab) ^ 0.65 day, and [Fe/H] < -1.6). NGC 6388 and NGC 
6441 (solid diamonds) have (Pab) values of 0.676 (nine a6-types) and 0.759 (42 a6-types), 
respectively . Note that NGC 6388 and NGC 6441 belong to neither Oosterhoff group. It is 
Pritzl et all J2000h who first showed that NGC 6388 and NGC 6441 do not fall into either 
of the usual Oosterhoff groups, and proposed the possibility of a new Oosterhoff class. Do 
NGC 6388 and NGC 6441 break down the traditional Oosterhoff classification scheme? Or, 
do we need to add a third class to the scheme? 



The answers to these questions depend solely on our understanding of the origin of the 
Oosterh off dichotomy. An exp lanation for t he ph enomenon has been put forward in the 



work of iLee fc ZinnI (119901 ) and lYoon fc Led (120021 ). In particular, lYoon fc Led (120021 ) have 
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discovered that most of the lowest-metalhcity ([Fe/H] < —2.0) GCs display a striking planar 
alignment in the outer-halo of the Milky Way. The alignment, combined with evidence from 
kinematics and stellar population, indicates that the metal-poorest GCs were originated from 
a satellite system, very likely from the Large Magellanic Cloud (LMC). In such a case, age 
and metallicity can be decoupled, so the lowest-metallicity Galactic GCs are not necessarily 
the oldest component of the Milky Way. This is interesting because the Oosterhoff dichotomy 
among the Galactic GCs can be naturally reproduced by assuming that the lowest-metallicity 
GCs are slightly younger than the genuine Galactic GCs of similar metallicity. 



Figure 56 conveys the essence of the lYoon &: Led (|2002| ) explanation. The Galactic GCs 
are re-classified according to their Galactocentric distances [Rg), and shown along with the 
population models. From this figure, one can see that the inner- {Rg < 8 kpc) and outer- 
halo {Rg > 8 kpc) GCs are well represented by the models of the 13-Gyr population (solid 
line), and by those of a slightly younger population (11.8 Gyr; dashed hne), respectively. 
The observed {Pab) of the lowest-metallicity ([Fe/H] < —2.0) GCs are well reproduced by 
the dashed line. This line also fits the GCs on retrograde orbits, whic h are most likely an 
accreted component of the Milky Way (Zinn 19931 : van den Bergh 1993). N o te als o that the 
LMC GCs follow the same dashed line. From this analysis. IYoou fc Lee ( 2002 ) conclude 
that the Oosterhoff dichotomy can be viewed as a manifestation of the age structure of the 
Galactic GC system, and thus, of the Galactic assembly history. 



Following the lYoon fc Led (|2002| ) framework, we now compare the AMR and SHR mod- 
els for NGC 6388 and NGC 6441 in terms of the {Pab) vs. [Fe/H] correlation (Fig. 56). The 
use of the input parameters in Table 3 leads the AMR model to predict an (Pgh) th at is too 



low at their [Fe/H] (the l ower stripe in Fig. 56). This result is in line with lPritzl et al.l (120021 ) 
and lCatelan et al.l (120031 ) . who confirmed the difficulty in reproducing the RR Lyrae periods 
for the given metallicities of NGC 6388 and NGC 6441. We have explored a large parameter 
space to find a way to reproduce the observations of NGC 6388 and NGC 6441 based on the 
AMR assumption, but failed unless we assumed an unreasonably high age (> 17 Gyr). Such 
an unrealistically old age is not acceptable. Moreover, the very old population is predicted 
to have too few RR Lyrae stars when compared with the observations. 

By contrast, the SHR model is remarkably consistent with the observations displayed 
in the Oosterhoff diagram (the upper stripe in Fig. 56) . We recall that all of the RR Lyraes 
belong to the blue-HB population in both the AMR and SHR models. When assuming that 
the helium abundance of the RR Lyraes is enhanced by AY = 0.05 (Table 3), their large 
(Pab) values can be easily achieved at their [Fe/H], even with RR Lyraes near the ZAHB. 
This is not surprising, since any effect that leads to a brighter HB suc h as an enhanced Y, 
would lead to longer periods, more consistent with the observations (ISweigart &: Catelan 
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19981 : ISweigartlll999f ) 



We finish this section by discussing the unusual values of N{c)/N{ab + c) of NGC 
6388 and NGC 6441 for t heir metaUicities. T he combined data of NGC 6388 and NGC 
6441 give 0.40 (= 34/85) (jCorwin et al.ll2006l ). The value falls between Oosterhoff group 
I (~ 0.2) and group II (~ 0.5). The c-type fraction of RR Lyrae variables has lon g been 
known to vary among GCs with differing HB types (ILee. Demarque. fc ZinnI Il990l ). The 
c-type fraction is significantly higher for predominantly blue HBs with little ZAHB stars 
wit hin the RR Lyrae strip ( Ooste rhoff group II) than for red or intermediate HBs (group 
I). iLee. Demarque. fc ZinnI (jl990l ) demonstrated that the evolutionary effect of blue HBs 
causes an uneven distribution in [B — V), since the speed of evolution increases as a star 
evolves further from its original position. As a result, the GCs with predominantly blue HB 
components are predicted to have high incidences of bluer RR Lyraes, i.e., c-type variables. 
In this regard, the GCs with 0.3 < N{c)/N{ab + c) < 0.4 can only be produced by a 
transitional HB morphology from intermediate- HB to predominantly-blue HBs. Under the 
normal-Y circumstances, GCs with such a transient HB type are rare, because it is hard for 
blue HBs to contain ZAHBs long enough to extend into the RR Lyrae strip. Few examples 
include NGC 5904 and NGC 6626 with A^(c)/iV(a6 + c) ~ 0.3. By constrast, the SHR model 
can generate relatively easily blue HBs with long ZAHBs extending into the instability strip. 
This is because the high-Y HB tracks are characterized by considerably stretched blue loops 
jSweigartl[l987l ). The unusual value of A^(c)/A^(a6 + c) (~ 0.4) for NGC 6388 and NGC 6441, 
therefore, points toward the SHR scenario. 

By means of the Oosterhoff class argument, we conclude that the SHR hypothesis pro- 
vides a more plausible solution to the problem posed by the unusual beh aviors of the RR 
Lyraes in NGC 6388 and NGC 6441. As suggested bv IPritzl et aP J200oh . NGC 6388 and 
NGC 6441 may be taken as a third class of Oosterhoff groups, and the helium abundance of 
RR Lyrae stars is the main driver behind this new class. 

In this section, we have compared the two models in terms of RR Lyrae properties 
(i.e., the pulsating periods and the Oosterhoff classes), and consistently found that the 
observations are strongly in favor of the SHR scenario over the AMR scenario. 



5. DISCUSSION 

Using stellar population simulations, we have investigated the possibility that the bi- 
modality of the HBs in NGC 6388 and NGC 6441 can be attributed to the presence of two 
distinct populations. This study was motivated by the discovery of multiple stellar popula- 
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tions coexisting in the most massive Galactic GC, uj Gen ( iLee et al.lll999bl : IStandford et al. 
20061 . and references therein). We are particularly ins pired by the fact that the integrated 
luminosities of NGG 6388, NGG 6441, and NGG 2 808 flHarrislll996h make them the 3rd, 5th , 
and 7th brightest among the ~ 150 Galactic GGs (lYoon. Lee, fc Ledl2000l : iRee et al.ll2002l ). 



Two chemical enrichment scenarios for GGs are discussed; namely, the age-metallicity 
relation (AMR) and the super-helium-rich (SHR) scenarios. In both models, two hypothet- 
ical populations coexist within individual GGs. We have examined the GMD morphologies 
and the RR Lyrae properties of NGG 6388 and NGG 6441, and found consistently that the 
SHR scenario is superior to the AMR hypothesis in accounting for the observations. Our 
best solution is that the blue HB and the RR Lyrae stars are comprised of helium-enhanced 
(AY ~ 0.05) stars, which make up ~ 15 % of the total population. Table 4 summarizes the 
model's ability to reproduce various aspects of the observations. 



Recently, iGaloi fc D'Antonal (120071 . hereafter GD07) have investigated the NGG 6441 
HB population, including the periods of the RR Lyrae variables. From their analysis, GD07 
conclude that Y ~ 0.35 is required for the tilted red clump, Y ^ 0.37 for the RR Lyrae stars, 
and Y ~ 0.38 - 0.40 for the blue HB. A partial, qualitative agreement is found between 
our SHR model and GD07's, in that the SHR population is necessarily invoked. However, 
several factors differentiate our result from GD07's: (a) While our model assumes two discrete 
populations with distinct helium abundances, GD07 present a continuous distribution of 
helium content; {b) The helium abundance required in our SHR model to reproduce the blue 
HB is Y ~ 0.29 - 0.30, which is well below the inferred Y in GD07 (Y ~ 0.35 - 0.40); (c) 
Our model exhibits the tilted red clump naturally reproduced by the inclusion of differential 
reddening effect, thus indicating the red clumps of a normal helium abundance (Y ~ 0.24 
- 0.25). In short, the main difference lies in the predicted total helium abundance of the 
GGs, in that our SHR model requires significantly lower value of total helium abundance. 
GD07 predict the total helium-enriched population amounting to ~ 60 % of th e entire stellar 
content. We note that, as mentioned in §§3.3, the R-method argument by iLayden et al. 



(119991 ) allows only a limited amount of helium cori t ent in the GGs. Furthermore, using 
deep HST photometry of NGG 6388, iGatelan et al.l (120061 ) discover the lack of a sizeable 
luminosity difference between the red HBs of NGG 6388 and 47 Tuc (a "flat" -red-clump 
GG with similar metallicity). This result points to the normal helium abundance of the red 
clump in NGG 6388, more consistent with our conclusion than GD07's. 

The high-heliu r n hyp othesis was originally suggested to explain the presence of hot EHBs 
by D'Antona et al. ( 2002h. This notion has been expanded to account for the bimodal HBs 



found in NGG 2808 (D'Antona fc Galoi 



(iD'Antona fc Galoill2004J : iD'Antona et al 



2004 



Lee et al 



2005 



2006 



2005h. NG G 6388, and NGG 6441 



Gatelanll2005l . GD07, this study). It is 
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noteworthy that the GCs with bimodal HBs often, if not always, also possess EHBs. Since 
the EHB phenomenon appears to share a common origin with HE bimodality, it is important 
to see whether the occurrence of EHBs correlates with GC mass, as does the instance of the 



bimodal HE (see Figure 1 of lLee. Gim. &: Casetti-Dinescul (120071 )). In fact, this was the basis 
for suggestions that the integrat ed luminosity plays an important role in the context of the 
second-parameter phenomenon (IFusi Pecci et al.lll993l ). Table 1 shows that eight Galactic 
GCs out of the 10 b rightest contain hot EHB components (the 7th column). According to 
Catelan et all JlQQsh catalog, ~ 90 % (14 out of 16) of the EHB GCs have blue HBs. Given 
the fact that increasing metallicity makes an HB redder, we speculate that the EHB GCs 
with normal blue HBs in Table 1 (NGC 6266, NGC 7078, and NGC 6273) could have typical 
bimodal HBs if they were as metal-rich as NGC 6388 and NGC 6441. 



Several mechanisms have been proposed to explain the peculiar chemical history of 



GCs with helium-excess populations ( 


3'Antona et al. 


2002 


: IVentu 


ra et al. 


2001, 


2002; 


Norris 


2004 




Gratton, Sneden, & Carretta 


2004 




3ekki & Norris 


20061: 


C 


luzhov 


200f 




Choi & Yi 


2007 


Tsujimoto, Shieevama, & Suda 


2007 




Newsham & Terndrup 


2007) 


. One theoretical 



requirement for the proposed chemical evolution is an initiall y high mass to retain the eiecta 



of the first gene r ation stars, as is evide nt in the case of u Gen ( lNorrisll2004j : iPiotto et al.ll2005 



NGC 6388, NGC 6441 flD'Antona fc Caloi 



Lee et al.l l2005l : iBekki &: Norrid 120061). The HB bimodality fouiid in raassiye GCs such as 



2004 



D'Antona et al 
2OO4I : 



2005 



Lee et al. 



2006; Catelan 



2005 



2005 



Piotto et al. 



CD07, this study), and NGC 2808 (iD'Antona fc Caloil 
20071 ) appears to be another piece of evidence for a link between the GC mass and the 
helium-enrichment process. It has been suggested that o o Gen was once part of a more 



massive system that later merged with the Milky Wav (ILee et al. 



Altmann. Catelan. fc Zoccalil l2005l : iMeza et al.l l2005l : iBekki fc Norrid |2006| . and references 



19993 : iDinescul 120021 



therein). Similar accretion events may have continued throughout the Galactic history, and 
the massive GCs with bimodal HB distributions may b e considered minor versions of a ; Gen , 
representing relics of the Galaxy assembly process (see iLee. Gim. fc Casetti-Dinescul (120071 ) 
for further discussion). 

Another important implication of our results concerns the long-standing second param- 
eter debate regarding HB morphology. Given that the HB-bimodality phenomenon has long 
been taken as the strongest evidence against the GC age being the second parameter, a 
fully satisfactory solution will not be found until the origin of the HB bimodality has been 
accounted for. Our results suggest that the GC-to-GC variation in heliur n abundance, which 
was one of the first candida tes for the second parameter in the literature (ISandage fc Wildey 
19671 : Ivan den Berghl Il967l ). have a role in giving rise to the HB bimodality. If our inter- 
pretation is confirmed, the GC-to-GC helium variation found most often in massive GCs 
should be considered a local effect rather than a global one, further supporting age as the 
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global second parameter. Further observations and modelling of GCs with bimodal HBs are 
still required to verify that helium variation is a local effect, and thus the third parameter 
controlling HB morphology. 

We would like to thank Michael Rich for providing the observational data in Figures 2 
& 3. Helpful comments from the referee were gratefully appreciated. S.-J. Y. acknowledges 
support from the Basic Research Program of the Korea Science & Engineering Founda- 
tion (Grant No. ROl-2006-000- 10716-0), and from the Korea Research Foundation Grant 
(MOEHRD: KRF-2006-331-C00134). 

REFERENCES 

Altmann, M., Catelan, M., & Zoccah, M. 2005, A&A, 439, 5 

Anderson, J. 2002, in u Centauri, a Unique Window into Astrophysics, ASP Conf. Ser. 265, 
eds F. van Leeuwen, J. D. Hughes, and G. Piotto. (San Francisco: ASP), p. 87 

Armandroff, T. E., & Zinn, R. 1988, AJ, 96, 92 

Beaulieu, S. F., et al. 2001, AJ, 121, 2618 

Bedin, L. R., Piotto, G., Anderson, J., Cassisi, S., King, I. R., Momany, Y., & Carraro, G. 
2004, ApJ, 605, L125 

Bekki, K., & Norris J. E. 2006, ApJ, 637, L109 

Bellazzini, M., Ferraro, F. R., & Buonanno, R. 1999, MNRAS, 307, 619 

Bellazzini, M., Fusi Pecci, F., Ferraro, F. R., Galleti, S., Catelan, M., & Landsman, W. B. 
2001, AJ, 122, 2569 

Brown, T. M., Sweigart, A. V., Lanz, T., Landsman, W. B., & Hubeny, I. 2001, ApJ, 562, 
368 

Busso, G., Piotto, G., & Cassisi, S. 2004, Mem. Soc. Astron. Itahana, 75, 46 
Caloi, v., & D'Antona, F. 2007, A&AS, 463, 949 (CD07) 
Castellani, V., & Quarta, M. L. 1987, A&AS, 71, 1 
Catelan, M. 2000, ApJ, 531, 826 



-17- 



Catelan, M. 2005, in Resolved Stellar Populations, ASP Conf. Ser., ed. D. Valls-Gabaud & 
M. Chavez, in press ( |asto-pli/0507464p 

Catelan, M., et al. 2006, ApJ, 651, L133 

Catelan, M., Borissova, J., Sweigart, A. V., & Spassova, N. 1998, ApJ, 494, 265 
Catelan, M., & de Freitas P. 1996, PASP, 108, 166 

Catelan, M., Sweigart, A. V., Pritzl, B. J., & Smith, H. A. 2003, in New Horizons in Globular 
Cluster Astronomy, ASP Conf. Ser. 296, eds. Giampaolo Piotto, Georges Meylan, S. 
George Djorgovski and Marco Riello (San Francisco: ASP), p. 289 

Choi, E., & Yi, S. K. 2007, MNRAS, 375, LI 

Chuzhoy, L. 2006, MNRAS, 369, L52 

Clement, C, et al. 2001, AJ, 122, 2587 

Clementini, G., et al. 2005, ApJ, 630, L145 

Corwin, T. M., Sumerel, A. N., Pritzl, B. J., Smith, H. A., Catelan, M., Sweigart, A. V., & 
Stetson, P. B. 2006, AJ, 132, 1014 

D'Antona, F., Bellazzini, M., Caloi, V., Pecci, F. F., Galleti, S., & Rood, R. T. 2005, ApJ, 
631, 868 

D'Antona, F., & Caloi, V. 2004, ApJ, 611, 871 

D'Antona, F., Caloi, V., Montaalban, J., Ventura, P., & Gratton, R. 2002, A&A, 395, 69 

D'Antona, F., Ventura, P., & Caloi, V. 2006, Mem. Soc. Astr. It., 77, 168 

Dinescu, D. I. 2002, in uj Centauri, a Unique Window into Astrophysics, ASP Conf. Ser. 265, 
eds F. van Leeuwen, J. D. Hughes, and G. Piotto. (San Francisco: ASP), p. 365 

Fusi Pecci, F., Ferraro, F. R., Bellazzini, M., Djorgovski, S., Piotto, G., & Buonanno, R. 
1993, AJ, 105, 1145 

Gratton, R. G., et al. 2007, A&A, 464, 953 

Gratton, R. G., Sneden, C, & Carretta, E. 2004, ARA&A, 42, 385 
Harris, W. E. 1996, AJ, 112, 1487 



Heitsch, F., & Richtler, T. 1999, A&A, 347, 455 
Hilker, M., & Richtler, T. 2000, A&A, 362, 895 

Kim, Y.-C, Demarque, R, Yi, S. K., & Alexander, D. R. 2002, ApJS, 143, 499 

Law, D. R., Majewski, S. R., Skrutskie, M. F., Carpenter, J. M., & Ayub, H. F. 2003, AJ, 
126, 1871 

Layden, A. C, & Sarajedini, A. 2000, AJ, 119, 1760 

Layden, A. C, Ritter, L. A., Welch, D. L., & Webb, T. M. A. 1999, AJ, 117, 1313 

Lee, H.-c, Yoon, S.-J., & Lee, Y.-W. 2000, AJ, 120, 998 

Lee, Y.-W., et al. 2005, ApJ, 621, L57 

Lee, Y.-W., Demarque, R, k Zinn, R. 1990, ApJ, 350, 155 

Lee, Y.-W., Demarque, P., & Zinn, R. 1994, ApJ, 423, 248 

Lee, Y.-W., Gim, H. B., & Casetti-Dinescu, D. I. 2007, ApJ, 661, L49 

Lee, Y.-W., Joo, J.-M., Sohn, Y.-J., Rey, S.-C, Lee, H.-C, & Walker, A. R. 1999b, Nature, 
402, 55 

Lee, Y.-W., Yoon, S.-J., Lee, H.-C, & Woo, J.-H. 1999a, in Spectrophotometric Dating of 
Stars and Galaxies, ASP Conf. Ser. 192, eds. Ivan Hubeny, Sally Heap, and Robert 
Cornett (San Francisco, ASP), p. 185 

Lee, Y.-W., & Zinn, R. 1990, in Confrontation between Stellar Pulsation and Evolution, 
ASP Conf. Ser. 11 (San Francisco, ASP), p. 26 

Lanz, T., Brown, T. M., Sweigart, A., Hubeny, I., & Landsman, W. B. 2004, ApJ, 602, 342 

Lejeune, T., Cuisinier, F., & Buser, R. 1998, A&AS, 130, 65 

Meza, A., Navarro, J. F., Abadi, M. G., & Steinmetz, M. 2005, MNRAS, 359, 93 

Moehler, S., & Sweigart, A. V. 2006, A&A, 455, 943 

Moehler, S., Sweigart, A. V., Landsman, W. B., Hammer, N. J., & Dreizler, S. 2004, A&A, 
415, 313 

Newsham, G., & Terndrup, D. M. 2007, ApJ, 664, 332 



-19- 



Norris, J. E. 2004, ApJ, 612, L25 
Oosterhoff, P. T. 1939, Observatory, 62, 104 

Pancino, E., Ferraro, F. R., Bellazzini, M., Piotto, G., & Zoccali, M. 2000, ApJ, 534, L83 
Park, J.-H, & Lcc, Y.-W. 1997, ApJ, 476, 28 

Piotto, G., ct al. 1997, in Advances in Stellar Evolution, eds. R. T. Rood and A. Renzini 
(Gambridge: Gambridge Univ. Press), p. 84 

Piotto, G., et al. 2005, ApJ, 621, 777 

Piotto, G., Bedin, L. R., Anderson, J., King, I. R., Gassisi, S., Milone, A. P., Villanova, S., 
Pietrinferni, A., & Renzini, A. 2007, ApJ, 661, L53 

Piotto, G., King, I. R., Djorgovski, S. G., Sosin, G., Zoccali, M., Saviane, I, De Angeli, F., 
Riello, M., Redo Blanco, A., Rich, R. M., Meylan, G., & Renzini, A. 2002, A&A, 
391, 945 

Pritzl, B., Smith, H. A., Gatelan, M., & Sweigart, A. V. 2000, ApJ, 530, L41 

. 2001, AJ, 122, 2600 

. 2002, AJ, 124, 949 

Pritzl, B., Smith, H. A., Stetson, P. B., Gatelan, M., Sweigart, A. V., Layden, A. G., & Rich, 
R. M. 2003, AJ, 126, 1381 

Raimondo, G., Gastellani, V., Gassisi, S., Brocato, E., & Piotto, G. 2002, ApJ, 569, 975 

Recio-Blanco, A., Aparicio, A., Piotto, G., De Angeh, F., & Djorgovski, S. G. 2006, A&A, 
452, 875 

Ree, G. H., Yoon, S.-J., Rcy, S.-G., & Lee, Y.-W. 2002, in Omega Gentauri: A Unique 
Window into Astrophysics, ASP Gonf. Ser. 265, ed. F. van Leeuwen, G. Piotto, and 
J. Hughes (San Francisco, ASP), p. 101 

Reimers, D. 1975, in Mem. Soc. R. Sci. Liege 6 Ser., 8, 369 

Rey, S.-G., Yoon, S.-J., Lee, Y.-W., Ghaboyer, B., & Sarajedini, A. 2001, AJ, 122, 3219 
Rich, R. M., et al. 1997, ApJ, 484, L25 
Rood, R. T. 1970, ApJ, 161, 145 



-20- 

Rood, R. T., et al. 1993, in The Globular Clusters-Galaxy Connection, ASP Conf. Ser. 48, 
eds. G. H. Smith and J. R Brodie (San Francisco, ASP), p. 218 

Rosenberg, A., Recio-Blanco, A, & Garcia-Marm, M. 2004, ApJ, 603, 135 

Sandage, A., & Wildey, R. 1967, ApJ, 150, 469 

Sarajedini, A., Chaboyer, B., & Demarque, P. 1997, PASP, 109, 1321 
Siegel, M. H., et al. 2007, ApJ, 667, L57 

Sosin, C, Dorman, B., Djorgovski, S. G., Piotto, G., Rich, R. M., King, I. R., Liebert, J., 
Phinney, E. S., & Renzini, A. 1997, ApJ, 480, L35 

Stanford, L. M., Da Costa, G. S., Norris, J. E., & Cannon, R. D. 2006, ApJ, 647, 1075 

Stetson, P. B., Vandenberg, D. A., & Bolte, M. 1996, PASP, 108, 560 

Sweigart, A. V. 1987, ApJS, 65, 95 

Sweigart, A. V. 1999, in Spectrophotometric Dating of Stars and Galaxies, ASP Conf. Ser., 
Vol. 192, ed. I. Hubeny, S. R. Heap, & R. H. Cornett (San Francisco: ASP), 239 

Sweigart, A. V., & Catelan, M. 1998, ApJ, 501, L63 

Sweigart, A. V., & Gross, P. G. 1976, ApJS, 32, 367 

Tsujimoto, T., Shigeyama, T., & Suda, T. 2007, ApJ, 654, L139 

van den Bergh, S. 1967, AJ, 72, 70 

van den Bergh, S. 1993, AJ, 105, 971 

van den Bergh, S. 1996, ApJ, 471, L31 

Ventura, P., D'Antona, F., & Mazzitelh, I. 2002, A&A, 393, 215 

Ventura, P., D'Antona, F., Mazzitelh, I., & Gratton, R. 2001, ApJ, 550, L65 

Yoon, S.-J., Lee, C.-H., & Lee, Y.-W. 2000, BAAS, 196, 41.02 

Yoon, S.-J., & Lee, Y.-W. 2002, Science, 297, 578 

Yoon, S.-J., Yi, S. K., & Lee, Y.-W. 2006, Science, 311, 1129 



- 21 - 



Zinn, R. 1993, in The Globular Clusters-Galaxy Connection, ASP Conf. Ser. 48, eds. G. H. 
Smith and J. P. Brodie (San Francisco, ASP), p. 38 



This preprint was prepared with the A AS IM^jX macros v5.2. 



-22- 



Table 1. The 10 Brightest Galactic Globular Clusters and the Features in their CMDs 



Rank 


NGC 


Name 




[Fe/H]'' 


Composite CMD<= ? 


HB Shape 




References 


1 


5139 


u) Ceil 


-10.29 


-1.59 


Yes (Multiple MS/HB/RGB) 


Multiple HB + 


EHB 


1, 2, 3, 4 


2 


6715 


M54 


-10.01 


-1.43 


Yes (Multiple RGB/HB) 


Multiple HB + 


EHB 


5, 6, 7 


3 


6388 




-9.82 


-0.60 


Yes (Bimodal HB) 


Bimodal HB + 


EHB 


8, 9, 10 


4 


2419 




-9.58 


-2.10 


No 


Normal blue HB 




8 


5 


6441 




-9.47 


-0.53 


Yes (Bimodal HB) 


Bimodal HB + 


EHB 


8, 9, 11 


6 


104 


47 Tuc 


-9.42 


-0.71 


No 


Normal red HB 




9 


7 


2808 




-9.36 


-1.37 


Yes (Multiple MS/HB) 


Bimodal HB + 


EHB 


9, 12, 13 


8 


6266 


M62 


-9.19 


-1.29 


No 


Long blue HB + 


EHB 


9 


9 


7078 


M15 


-9.17 


-2.17 


No 


Normal blue HB 


+ EHB 


9 


10 


6273 


M19 


-9.08 


-1.68 


No 


Normal blue HB 


+ EHB 


9 



''M*?* is obtained from lHarrij l ll996h . 
''[Fe/H] is obtained from lZinnI lll993l '). 
'^CMDs with multiple RGBs and/or HBs. 



References. — (1) iLee et al.l lll99 9tj): (2) 'Hilker & Richtlei' ('2OOO'); (3) 'Pancino et al.l l!2000l) ; (4) JPiotto et a ll l'2005');_(5) 
iBcllazzini ct al.llll999h: f el lLavden fc Saraic dini (2000'): (7) Rosenberg et al. (20 04): (81 Rich et al. (Il997l ): (Q'l lPiotto et al., (.2002) : 
aO) ICatelan et al.l (I2OO6I) : f 11) IPritzl et al.l J2OO3I) : (12) ISosin et aLI (ir997h : ('13) IPiotto et al.l (l2007l) 
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Table 2. Model Ingredients 



Ingredient Stellar & Flux Libraries 

MS to RGB Evolutionary Tracks Kim et al. (2002) Yonsei-Yale Normal-helium Isochrones 

Kim et al. (2008, in prep.) Yonsei-Yale Helium-rich Isochrones 

Post-RGB Evolutionary Tracks . Han et al. (2008, in prep.) Yonsei-Yale Normal-helium HB Tracks 

Han et al. (2008, in prep.) Yonsei-Yale Helium-rich HB Tracks 

Flux Library Leieune. Cuisinier. fc Buser (1998) Model Atmosphere 

RR Lyrae Variable Stars Lee. Demarque. fc Zinn (1990) Prescriptions 
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Table 3. Model Input Parameters 



Parameter 



NGC 6388 



NGC 6441 



Common Parameters: 

Initial mass function Salpeter 

a-element enhancement, [a/Fe] 0.3 

Distance modulus, {V — My) (mag.) 16.50 

Galactic reddening, E{B - V) (mag.) 0.34 

Differential reddening, crE(B-v) (mag-) 0.03 

Extinction coefficient, Ay /i^(i3 — y) 3.1 

HB mass dispersion, ctm (Mq) 0.02 

RR Lyrac instability strip width, ALogTe (K) 0.085 

Rcimcrs' (1975) mass-loss efficiency parameter, rj 0.56 

The AMR ModeP : 

Metal abundance, Z 0.007 & 0.0035 

Helium abundance, Y 0.245 & 0.245 

Absolute age, t (Gyr) 13.0 & 15.0 

Number fraction (%) 82 & 18 

The SHR Model^^ : 

Metal abundance, Z 0.007 & 0.007 

Helium abundance, Y 0.245 & 0.295 

Absolute age, t (Gyr) 13.0 & 13.0 

Number fraction (%) 82 & 18 



Salpeter 

0.3 
17.17 

0.43 

0.03 

3.1 

0.02 
0.085 

0.56 



0.008 & 0.0055 
0.245 & 0.245 
13.0 & 15.0 
86 & 14 



0.008 & 0.008 
0.245 & 0.295 
13.0 & 13.0 
86 & 14 



^Metal-rich younger & metal-poor older populations 
'^Normal-helium & helium-enriched populations 
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Table 4. Summary of Model Ability to Reproduce Observations 



Feature 


The AMR Model 


The SHR Model 


Single RGB 


O 


o 


RGB Bump Luminosity 


o 


o 


RGB Bump Slope 





o 


HB Bimodality 


o 





Red HB Tilt 








Overall HB Slope 


X 





Extreme HB 


X 


A 


RR Lyrae Luminosity . . 


X 


o 


RR Lyrae Period 


X 


o 


N{c)/N{ah + c) 


o 


o 


Oosterhoff Diagram 


X 


o 
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Fig. 1 . — (a) The V^-band luminosity function (open histogram) of 147 Galactic GCs from 



Harris! (119961 ) catalog. Arrows denote the integrated magnitude of the brightest GC, u Gen, 



and the GGs with bimodal HBs among the 10 brightest GGs. The numbers in the parentheses 
denote their brightness ranks. NGG 6388 and NGG 6 441 are among themost massive GGs. 
The filled histogram is for the 11 bimodal-HB GGs (IGatelan et al.l Il998l ) plus u Gen and 
M54. Note that they are all brighter than My* ^ —7. (6) The relative number fraction 
of bimodal-HB GGs plus u Gen and M54 as a function of the integrated magnitude. Their 
fraction in a luminosity bin increases with increasing luminosity. Error bars show Poisson 
errors. 
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Fig. 2. — Comparison of the observational data (a & d) with the AMR model {b & e) and 
the SHR model ( c & /) for NGC 6388 and NGC 6441. (a & d) The observed CMDs are from 



Rich et al.l (119971 ). Isochrones are overlaid with observations, from which the dominant major 
populations, i.e., the metal-rich younger population in the AMR model and the normal- 
helium populatio n from the SHR naodel, are ge nerated. Crosses represent the observed RR 
Lyrae stars from iPritzl et al.l (l200ll . 120021 . |2003| ). The same RR Lyrae stars are also shown 
in Fig. 4. Horizontal dashed lines represent the brightness at the base of red clumps. Short 
tilted solid lines denote the RR Lyrae instability strips. (6, c, e, & /) In 6 & e, the metal- 
rich younger (denoted by "1") and the metal-poor older (denoted by "2") populations are 
respectively shown as red and blue dots. The magenta and cyan isochrones are respectively 
for the metal-rich younger and the metal-poor older populations. The combinations of 
metallicity, helium abundance, and age required to reproduce the CMD morphology are 
denoted. Note that the observed bimodal HBs are explained as the sum of the red HBs from 
component 1 (red dots) and the blue HBs from component 2 (blue dots). Crosses denote the 
model RR Lyrae stars. In c & /, the normal-helium (denoted by "1") and the helium-excess 
(denoted by "2" ) populations are respectively shown as red and blue dots. The magenta and 
cyan isochrones are respectively for the normal-helium and the helium-excess populations. 
Note that the observed bimodal HBs are explained as the sum of red HBs from component 1 
(red dots) and blue HBs from component 2 (blue dots). Crosses denote the model RR Lyrae 
stars. 
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Fig. 3. — Numerical experiments on the effect of differential reddening on the red clumps. 
The observed clumps (a & d) are compared with the fiat clumps generated without the 
differential reddening effect (6 & e), and the tilted clumps generated with the effect (c & /). 
Note that panels (a & d) and (c & /) arc identical to Figs. 2a & 2d and Figs. 2c & 2/, but 
have been magnified around the HB red clump region. The dominant major population is 
common across the AMR (the metal-rich younger component denoted by "1" in Fig. 2) and 
SHR (the normal-helium component denoted again by "1" in Fig. 2) models. The inferred 
differential reddening appears to be enough to turn a "normal" red HB component into a 
significantly sloped structure. Note also that the differential reddening effect appears to be 
necessary to reproduce the observed scatter and the possible tilt of the RGB bumps. 
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Fig. 4. — Comparison of the observations (a & d) with the AMR model (6 & e) and the SHR 
model (c & /) for NGC 6388 and NGC 6441 in the (B-V {Pp)) plane, (a & The observed 
RR Lyrae variables are from lPritzl et al.l (l200ll . |2002| . 120031 ) . Solid and open squares represent 
ab-tjpe and c-type RR L yrae stars, respectively. Type c RR Lyrae were fundamen talized by 



Log Pp = Log Pc + 0.13 (ICastellani &: Quarta 



1987; 



Lee. Demarque. fc Zinrull990[). In order 



to avo id small number statistics, the catal ogs of NGC 6 388 (12 RR Lyraes from lPritzl et al. 



((20021)) and NGC 6441 (63 RR Lyraes from lPritzl et al.l (|200l|,l2003|)) were combined. Dashed 
line in (a) denotes the observed mean Pp value {{Pp) = 0.661 day) for the combined NGC 
6388 (magenta) and NGC 6441 (cyan) data. Solid line is the least-squares fit to the data. 
Panel (d) is identical to panel (a). (6, c, e, & /) Examples of the Monte Carlo simulations. 
Solid lines are the same as in {a Sz d). Dashed and dotted lines in {b & e) and (c & /) denote 
the mean values of Pp and their 1 a uncertainties of 10,000 synthetic simulations for the 
AMR model and the SHR model, respectively. 
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Fig. 5. — (a) Oosterhoff classes. The correlation between the mean period of type ab RR 



Lyrae s {{Pab}) and [Fe/H] for GCs in the Galaxy is shown. The data are from I Clement et al. 



(120011 ). GCs having 5 or more type ab RR Lyraes are shown. The error bar denotes the 
mean errors. GCs are mainly divided into two distinct groups according to the mean period 
and the metal abundance: Oosterhoff group I (filled squares; (Pab) — 0.55 days, and [Fe/H] 
> —1.8) or group II (open squares; {Pab) — 0.65 days, and [Fe/H] < —1.6). The {Pab) 
values for NGC 6388 (0.676 days from nine typ e-afe variables) and fo r NGC 6441 (0.759 days 
from 42 type-a6 variables) are obtained from Corwin et al. ( 20061 ). Note that NGC 6388 
and NGC 6441 (solid diamonds) belong to neither Oosterhoff group. (6) Inner- & outer-halo 
division. Solid and dashed lines are our model predictions for the 13-Gyr and 11.8-Gyr 
populations, respectively. GCs with the Galactocentric distance, Rq < 8 kpc (filled squares) 
are well reproduced by the models for the old population (solid line), whereas most GCs 
with Rc > 8 kpc (open squares) follow the model locus for slightly younger ages (dashed 
line). Small plus signs mark GCs on retrograde orbits around the Galactic center, whereas 
triangles represent the old GCs found in the LMC halo. Note that the retrograding GCs 
and the LMC GCs follow the dashed line. Large gray cross between the diamonds marks 
the number-weighted mean values of [Fe/H] and {Pab) for NGC 6388 and NGC 6441. Two 
thick stripes indicate the AMR model (At = +2 Gyr, lower stripe) and the SHR model (AY 
= +0.05, upper stripe). Note that the SHR model prediction succeeds in reproducing the 
unusual positions of NGC 6388 and NGC 6441. 



